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ABSTRACT 

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) lumi¬ 
nosity function does not significantly evolve at 4 < z < 7 and is approximately ^-21. We investigate 
this apparent non-evolution by examining a sample of 178 bright, A/yy < -21 galaxies at z = 4 to 7, analyzing 
their stellar populations and host halo masses. Including deep Sp/fzer/IRAC imaging to constrain the rest-frame 
optical light, we find that Mj^y galaxies at z = 4-7 have similar stellar masses of \Qg{M/M q) = 9.6-9.9 and are 
thus relatively massive for these high redshifts. However, bright galaxies at z = 4-7 are less massive and have 
younger inferred ages than similarly bright galaxies at z = 2-3, even though the two populations have similar 
star formation rates and levels of dust attenuation for a fixed dust-attenuation curve. Matching the abundances 
of these bright z = 4-7 galaxies to halo mass functions from the Bolshoi ACDM simulation implies that the 
typical halo masses in ^ Myy galaxies decrease from log(Mh/M 0 ) = 11.9 at z = 4 to log(Mh/M 0 ) = 11.4 at 
z = 7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live 
in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Hb/ ftm rises 
from 5.1% atz = 4 to 11.7% atz = 7; this evolution is significant at the ^ 3a level. This rise does not agree 
with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency 
of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts. 

Subject headings: early universe — galaxies; evolution — galaxies; formation — galaxies; high-redshift — 
ultraviolet; galaxies 


1. INTRODUCTION 

Tracing the buildup of stellar mass from the epoch of the 
first galaxies through the present can be used to constrain 
models of galaxy formation. Understanding what physics 
governs this buildup is one of the key outstanding questions 
in galaxy evolution. There is a consensus that the cosmic star 
formation rate (SFR) density rises from the dawn of galaxies, 
peak s at redshifts z 2-3. and t hen declines steeply at z < 2 
(e.g.. lMadau & Dickinsonll2014h . The physical origin of this 
evolution in galaxy stellar mass growth is poorly understood, 
especially at high redshifts. While studies of galaxy evolu¬ 
tion routinely quantify the stellar content of distant galaxies, 
it remains challenging to relate the stellar masses to the sup¬ 
ply of gas fueling star formation. Theoretical works attempt 
to address these fundamental questions to a varying degree 
of success, but observational data have remained incomplete, 
particularly at the massive end, and they are dominated by 
systematic uncertainties unavoidable in the stitching together 
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of datasets from different campaigns (e.g., Behroozi et al. 
2013a). 

The advances facilitated by the Wide Field Camera 3 
(WFC3) on the Hubble Space Telescope (HST) over the past 
half-decade hav e led to the discovery of >100 0 redshift z > 
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facilitated by these large samples is the measurement of the 
rest-frame ultraviolet luminosity function, which quantifies 
the relative abundances of galaxies over a wide dynamic range 
in luminosity. As the UV light probes recent star formation 
activity, the integral of the rest-frame UV luminosity func- 
tion provides an estimate of the cosmic SFR density (e.g.. 


Madau et al.lll996tlBouwens et al.ll^l2HMadau & DickinsonI 


20141; iFinkelstein et al.l [2015h . The lumin osity function is 

typically parameterized with the ISchechte^ (^1976^ functional 
form that is a power law at low luminosities and declines ex¬ 
ponentially at high luminosities. Its parameters are the char¬ 
acteristic luminosity M^y, the faint-end slope a, and the nor¬ 
malization (jf. Previous studies typically found that these 
parameters evolved with redshift; the characteristic luminos¬ 
ity decreased w ith increasing redshift as the f aint-end slope 
steepened (e.g., iBouwens et al.ll2007l 1201 lallbl ; iMcLure et al.l 
120131) . This “luminosity evolution” of the luminosity function 
was widely accepted as it fit the general trend observable in 
the evolution of the cosmic SFR density. 

More recent work, however, has shown that the picture de- 
scribed above is i ncomplete. The first ev idenc e came from 
lOno et all (120121) . IFinkelstein et al.l (120131) . and iBowler et al.l 
(120141) . where a larger than expected number of bright galax¬ 
ies turned up as z = 7 surveys moved to wider fields. 
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Fig. 1.— The spectral energy distribution of a representative galaxy in each of our redshift bins, shown as the blue circles. For filters where the signal-to-noise 
was <2, we show the Icr upper limits. The black curve shows the best-fit stellar population model, and we list the ID and best-fit stellai' mass for each galaxy. 


iFinkelstein et alJ (120151) and iBouwens et all (120151) have con¬ 
firmed this excess, and recent studies h ave concluded it is 
not attributable to gravitati onal len sing (iMason et'd] 201^ 
iBarone-N ugent et 'ffl l2015h . Both IFinkelstein et alJ ( 201^ 
and iBouwens et all (120151) computed the evolution of the rest- 
frame UV luminosity function at z = 4-8, finding that con¬ 
trary to the preceding results derived from smaller datasets, 
the characteristic luminosity was remarkably redshift- 
independent between z = 4, 5, 6, and 7. The constancy broke 
down only at z = 8 where the data were least constraining. 
Most of the evolution took place in the characteristic number 
density: it declined towards higher redshifts. Therefore, while 
galaxies in general became less common at higher redshifts— 
consistent with the decline in the cosmic SFR density—^bright 
galaxies remained relatively common in the distant universe. 

Here, we seek to constrain the physical properties in dis¬ 
tant UV-bright galaxies and attempt to understand how they 
maintained high levels of star formation. In §2, we de- 
scribe the brigh t galax y sample that we have taken from 
IFinkelstein et al.l (120151) and discuss the additional constraints 
that can be placed on the stellar populations with the inclusion 
of Sp/fzer/IRAC photometry. In §3, we use cosmic abundance 
matching to estimate the halo masses and the stellar-to-halo 
mass ratios for these galaxies. In §4, we discuss the evolu¬ 
tion of the stellar baryon fraction with redshift, and in §5, 
we present our concl usions. We assume the WMAP7 ACDM 
cosmological model (iKomatsu et akllMITl) throughout, with 
Hq = IQ.2 km s-' Mpc-‘, = 0.275, and Ha = 0.725. A ll 
magnitudes given are in the AB system (lOke & Gunnlll98^ . 

2. STELLAR POPULATIONS IN UV-BRIGHT GALAXIES 

Here we wish to constrain the physical processes that reg¬ 
ulate the abundance of bright galaxies in the distant universe. 
These luminous systems are observed only a short time af¬ 
ter the Big Bang and trace prominent density peaks at their 
epoch. We further investigate sources withMisoo < -21. This 
magni tude is approximately the value of M*^y at th ese red¬ 
shifts (iBouwens et al.ll^l5t IFinkelstein etd] 120151) . though 
the exact value of does get progressively more uncer¬ 
tain wi th increasing redshift ( ifc 0.09 at z = 4 to ± 0.4 at 
z = 7: IF inkelstei n et al.l 120151) . Using ground-based data, 
iBowler et al.l (120151) have found evidence that M^y may be 
fainter than -21 at z = 7, though only at the 2tT level, and 
thus not significantly disc repant with the measur ements of 
IBouwens et al.l (120151) and IFinkelstein et al.l (120151) . Because 
this study is concerned with the physics driving the apparently 
high star-formation rates in distant bright galaxies, the exact 
luminosity we choose is not critical. In order to explore rel¬ 
atively bright galaxies, we choose the approximate value of 


Myy at these redshifts as our threshold. 

We use the sample o f ^ 7500 galaxies at 4 < z < 8 from 
IFinkelstein et al.l (120151) and refer the reader to that paper for 
details on the photometry, photometric redshift sample se¬ 
lection, and derivation of UV absolute magnitude at 1500 A 
(Muv)- From their full catalog, here we analyze the 150, 75, 
28, 18, and 3 galaxies with Muv < -21 at z = 4, 5, 6, 7, and 
8, respectively, which come from the CANDELS GOODS- 
South and North fields (with a total area of ^ 280 arcmin^). 
The redshift bins are bounded by Az = ±0.5 from the central 
redshift. 


2.1. Inclusion of Spitzer/IRAC Photometry 

To learn more about these intriguing systems, we turn to 
stella r population modeling. We use the HST photometry 
from IFinkelstein et al.l (120151) that includes ACS and WFC3 
PSF-mat ched total fluxes in the wavelength range 0.4 - 1.8 
/rm (see iKoekemoer et al.ll201 ll for details on the imaging). 
This catalog also includes photometry fro m the Spitzer Space 
Telescope Infrared Array Camera (IRAC; iFazio et ^I2004^ 
imaging of our fields. The IRAC imaging probes the rest- 
frame optical at z > 4 and thus provides significant con¬ 
straining power on the stellar masses of our galaxies. It also 
gives a more accurate handle on the ages and dust attenu¬ 
ations by reducing the degeneracy between the two param¬ 
eters. The details of the long-wav elength photometry are 
presented in IFinkelstein et al.l (120151) : here, we review them 
only briefly. The mosaics were obtained by coadding all the 
available data in these fields: the GOODS (Dic kinson et al., 
in pre p), Spitzer Extended Deep Survey (SEDS; Ashb^_et^ 


2013 ). and Spftzer-CANDELS (S-CANDELS; 


Ashby et al 


20151) wide-held programs, as well as the deep pointings 


from Spitzer prog ram 70145 (the IRAC Ultra-Deep Eield of 
iLabbeet al.ll2^ over the Hubble Ultra Deep Eield and its 
parallels, and program 70204 (PI Eazio) which observed a re¬ 
gion in the GOODS-S held to 100 hr depth. The hnal mosaics 
have a depth of > 50 hr over both CAN DELS GOODS hel ds 
and >100 hr over the H UDE main held (lAshbv et al.ll2015h . 

The TPHOT software (iMerlin et aLlUOlsI) . an updated ver¬ 
sion of TFIT (Laidler et al. 2007), was used to measure 
photometry in the Spitzer/IRAC imaging. This software 
models the low-resolution IRAC images by convolving the 
//577WEC3 //-band image with an empirically derived IRAC 
PSE, simultaneously htting all IRAC sources. This provides 
robust photometry even for moderately blended sources. The 
full descript ion of our TPHOT IRAC photometry catalog is 
presented in ISong et akl (120151) . 

All high-redshift galaxies were visually inspected in the 
TPHOT residual maps. If an object was on or near a strong 
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residual, reliable IRAC photometry was not possible, affect¬ 
ing 20-30% of the galaxies in our bright galaxy sample. To 
obtain the most robust stellar mass measurements, in our sub¬ 
sequent analysis we do not include these affected galaxies. 
Over 90% of the remaining galaxies in our sample had a 3.6 
pm or 4.5 pm detection of at least 3a significance, with a 
magnitude range at z > 6 of 23.5 < OT3.6 < 25.5. This is 
expected, as galaxies with M 1500 < -21 should be massive 
enough at all redshifts to yield an IRAC detection absent 
crowding. In fact, when comparing median stellar masses 
derived excluding and including galaxies without IRAC con¬ 
straints, the median stellar mass of galaxies with IRAC con¬ 
straints is at most 0.1 dex lower than that of the whole sam¬ 
ple. This is because galaxies with true mid-infrared fluxes 
well below the IRAC detection limit can have poorly con¬ 
strained SEDs. In contrast, our results show that the typical 
UV bright z > 6 galaxy has a lower mass-to-light ratio than 
other possible solutions, thus the IRAC detection prior does 
not drive us to higher M/L models. 

The IRAC photome try used here is the same as by 
iFinkelstein et al.l (120151) who used the IRAC and additional 
HST/ACS F814W photometry to re-measure the photometric 
redshifts, removing 14, 14, and 1 galaxies from their z = 4, 
5, and 6 samples, respectively. The galaxy sample we con¬ 
sider h ere comes from the cleaned sample of iFinkelstein et ahl 
(120151) . thus these presumed lower redshift interlopers have 
already been removed. In this work, we have performed an 
additional iteration of visual inspection of the IRAC imag¬ 
ing, which results in a few more galaxies having preferred 
lower redshift solutions. From our sample we remove seven 
such additional galaxies which now have preferred lower red¬ 
shift solutions (two from our z =4 sample, and five from our 
z =5 sample). We also remove five additional sources from 
our z = 4 sample which have 2.7 < Zphot < 3.5, as we do not 
wish to bias our z = 4 sample stellar mass measurements. Af¬ 
ter excluding these objects, galaxies in the remaining sample 
all have photometric redshifts within Zsampie - 0.5 < Zphot < 
^sample + 0.7, thus we Consider the effects of potential low- 
redshift interlopers to be minimal. The inclusion of the very 
small number of sources at Zphot = 2sampie+ 0.5 - 0.7 (eight at 
z = 4, six at z = 5 and 1 at z = 7; see tables in Appendix) make 
no difference to the med ian stellar mass discussed below. See 
IFinkelstein et al.l (120151) for a quantitative description of the 
potential contamination. 

2.2. Spectral Energy Distribution Fitting 

The technique we used to fit stellar population models 
to photometry was sim il ar to th e one we emp loyed before 
(IFinkelstein et all |2^ l2012bllaL IMl IMTl . We used 
the updated (2007) stella r population synthesis models of 
iBruzual & Charlo3 (1200.31) to generate a grid of model spec¬ 
tra*. We varied the stellar mass (defined as the total gas 
mass converted into stars), the stellar population metallicity, 
the time since the onset of star formation (henceforth, the 

* T he 2007 update to ttie stellar population models of lBruzual & Charloll 
1200.31) may overestimate the contribution of thermally pulsating asymptotic 
giant branch (TP-AGB) stars. However, these stars typically begin to dom¬ 
inate the emission at population ages >1 Gyr and rest frame wavelengths 
> 1 ^m. Our longest wavelength filter (4.5 /rm) at our lowest redshift {z = 
4) probes only 0.9 /rm, and all other filter/redshift combinations probe bluer 
rest-frame wavelengths. The TP-AGB contribution may impact the SEP i n 
post-starburst galaxies at wavelengths as low as 0.5 /rm iKriek et al.ll20101) . 
However, our galaxies are highly star-forming, with inferred population ages 
are <C 1 Gyr. Thus, our choice to use the updated models likely has no effect 
on our results. 


age), and the star formation history (SFH). We assumed the 
Salpeter^ initial mass function (IMF). Allowed metallicities 
spanned (0.02 - 1)Zq and ages spanned 1 Myr to the age 
of the universe at the source redshift. We allowed several 
different SFH scenarios, including a single burst, continuous 
star formation, and both the exponentially decaying and ris¬ 
ing (so-called “tau” and “inverted-tau”) models. We included 
nebula r emission lines using the prescri ption ofISalmon et all 
(120151) . which takes the line ratios from llnou3 (1201 ll) . assum- 
ing that the gas has the same metallicity as the stars and that 
all the ionizing photons emitted by the model stellar popu¬ 
lation are reprocessed in the galaxy and their escape is neg¬ 
ligible. To the rest frame spectra we a dded dust attenuation 
using the starburst attenuation curve of ICalzetti et alJ (120001) 
in the range of 0 < E(B-V) < 0.8 (0 < Ay < 3.2 mag). Then 
we redshifted the models to 0 < z < 11 and ad ded intergalac- 
tic medium (IGM) attenuation (lMadaulll995l) . The resulting 
model spectra were integrated through our HST and Spitzer 
Alter bandpasses to derive synthetic photometry for compari¬ 
son with our observations. 

We emphasize that our model parameterization assumes 
that the SFH of each object follows one of the simple sce¬ 
narios (i.e., single burst, continuous, tau, or inverted-tau), not 
a superposition of such scenarios. This may seem like an 
oversimplification, but evidence is mounting t hat the SFHs in 
distant galaxies vary smoothly with time (e.g.lPapovich et al. 
2011 : l^lator et all 1201 It iReddv et^l2012bt ISalmon et al. 
201-51) . Therefore, the simple scenarios may in fact be rather 
good appro ximations, in pa rticular when deriving the stellar 
mass (e.g.. iFee et akllMI^ . Therefore, the SFHs and thus 
the physical properties of the bright galaxies that we con- 
sider here should not be strongly affected by burstiness (e.g., 
iJaacks et al.ll2012l) . 

The best-fit model was found via minimization. We in¬ 
cluded an extra systematic error of 5% of the object flux in 
each band to crudely account for the residual uncertainties in 
the zero point correction and PSF matching process. The un¬ 
certainties in the best-fit parameters were derived via Monte 
Carlo simulations, perturbing the observed flux of each object 
in each Alter with a number drawn from a Gaussian distribu¬ 
tion with a standard deviation equal to the flux uncertainty in 
the Alter. Taking the source redshift to be statistically uncor¬ 
related with other spectral energy distribution (SED) fitting 
parameters, in each Monte-Carlo realization we drew the red¬ 
shift from the photometric redshift statistical likelihood func¬ 
tion of the given object. To prevent low-redshift solutions 
from biasing the physical parameters, we limited the random 
redshift to be within Az = ± 1 of the best-fit photometric red¬ 
shift. This treatment of the source redshift effectively folded 
the uncertainty in redshift into the uncertainty in the phys¬ 
ical parameters (most notably, the stellar mass and M 1500 ; 
IFinkelstein et al.ll2012al) . For each galaxy, 10^ Monte Carlo 
realizations were generated and this provided a sample of as 
many values for each model parameter. SFRs for the best-fit 
model and for each Monte Carlo realization, were derived by 
converting the dust-correcte d value of M 1500 to a SFR via the 
relation of iKennicutll ([l99^. 

For our subsequent analysis, we discarded realizations with 
poor best-fit models with > 20 to ensure robustness of the 
derived properties. This removed a relatively small number of 
galaxies (16, 9, 1,0 and 1 at z = 4, 5, 6 , 7 and 8 , respectively). 

^ To convert our results to tliose obtained from a Cliabrier IMF, one sliould 
divide the stellar masses and SFRs by a factor of 1.7. 
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Fig. 2.— Joint probability distributions for the stellar masses of bright 
galaxies (M 1500 < —21) in our z = 4, 5, 6 and 7 galaxy samples (blue line). The 
histograms are of the best-fit values, and the red dashed line and red shaded 
regions denote the median and Icr uncertainty on the median, respectively, of 
the joint PDFs. Little significant evolution in the median stellar mass is seen 
over this redshift range. 

The final sample contains 94, 46, 19, 14, and 1 galaxy at z = 
4, 5, 6 , 7, and 8 , respectively. As there is only one galaxy at 
z = 8 that satisfies these criteria, we focus on 4 < z < 7 for 
the remainder of this paper. Figure [T] shows the SED fit for a 
typical galaxy in each of our redshift bins. 

2.3. Physical Properties 

The results of SED fitting are summarized in Table [T] To 
derive the median stellar population properties, rather than 
stacking the images or fluxes of the galaxies, we stacked 
samples of Monte Carlo realizations in the model parameter 
space. In each redshift bin, the stacked sample allowed us 
to quantify the multivariate distribution of galaxy properties. 
The joint probability distribution functions for stellar mass 
marginalized over other parameters are shown in Eigure |2l 
The median of a parameter such as the stellar mass is taken 
to be the median value in the bin. The Icr confidence interval 
on the median was calculated via 10 ^ bootstrap simulations 
where we rederived the median from a randomly drawn (with 
replacement) sample of the galaxies in each redshift bin. 

As shown in Table [T] the median stellar population param¬ 
eters do not evolve significantly with redshift. Broadly speak¬ 
ing, Muv < -21 galaxies at z = 4-7 are moderately massive 
(log[M,/M 0 ] Ri 9.6-9.9), somewhat young (<100 Myr), and 
have non-negligible dust attenuation {E(B-V) = 0.07-0.13 
corresponding to Ay w 0.3-0.5 mag) and high SERs (^40-60 
Mq yr-i). 

Significant amounts of dus t are likely produced in g alax¬ 
ies as early a s z 7, as lEinkelstein et al.l (12012bl) and 
iBouwens et al.l (l2014h have previously noted that massive 
and/or UV-bright galaxies had similarly red UV continuum 
slopes at z = 4-7, with a typical value of the UV spectral slope 
/3 ^ -1.8. Both studies concluded that this implied a similar 
amount of dust in bright/massive z = 4-7 galaxies, indepen¬ 
dent of redshift. We confirm this result, finding £’(Z?-y) ~ 0.1 
in bright galaxies at z = 4-7 (constrained at 68 % confidence 
to be > 0 and < 0.15). Therefore, although fainter/1 ower- 
mass galaxies appear to be less dusty at highe r redshifts (e.g., 
lEinkelstein et aljf2012btlBouwens et alJl2014ft . this is not true 
for the brightest galaxies. This can be confirmed with ALMA, 
and in fact ALMA has recently detected dust emission from 
normal galaxies out to z 5-7.5 (e.g.. I Watson et al] 1201 51 : 


TABLE 1 

Median Physical Properties of Galaxies with M 1500 < -21 


Reds (lift 

Number 

log(M*/M 0 ) 

Age 

(Myr) 

E(B-V) 

SFR 

(Mq yr-‘) 

z = 4 

94 

9.86 ± 0.04 

44 ±2 

0.13 ±0.01 

56 ± 4 

z = 5 

46 

9.80 ± 0.06 

35 ±2 

0.12 ± 0.02 

52 ± 10 

z = 6 

19 

9.78 ± 0.07 

40 ±4 

0.07 ± 0.02 

40 ± 8 

z = 7 

14 

9.64 ±0.13 

29 ±8 

0.09 ± 0.02 

41 ±9 


Note. — Median values of physical parameters from the joint prob¬ 
ability distribution describing all galaxies in a given redshift bin that 
have a measurement in the IRAC 3.6 pm channel and a best-fit model 
with < 20. The number of galaxies satisfying these criteria in each 
redshift bin is given in the second column. The statistical uncertain¬ 
ties on median values were derived via 10^ Monte-Carlo simulations in 
which the median was rederived from a randomly drawn sample (with 
replacement) of the galaxies in each redshift bin. The spread in values 
for individual galaxies is larger. 

ICapak et'^l2015[) . 

If the amount of UV attenuation due to dust among bright 
galaxies had evolved from z = 7 to 4, then it could have led 
to our selecting a lower stellar mass at a given UV absolute 
magnitude (a similar effect could occur if the dust attenuation 
curve is redshift dependent). However, our results show that 
this is not the case, as not only does our inferred attenuation 
exhibit no evolution, but the median stellar mass also appears 
roughly constant over the redshift interval. 

Our fiducia l analy sis assumed a dust attenuation curve from 
iCalzetti et al.l (120001) for consistency with previous results in 
the literature. However, a number of recent studies ha ve found 
that a Small Magellanic Cloud (SMC)-like (iPeil 19921) attenua¬ 
tion curve may be more appropriate for high-redshift galaxie s 
(e.g.. l^ddv et al.l2012al:lTilvi et al.l2013llCapak et al.l20L5l) . 
This dust attenuation curve has more attenuation for given 
values of E{B-V), thus we explored how our results would 
change had we assumed this attenuation curve in our fidu¬ 
cial analysis. Even with this different attenuation curve, the 
amount of dust attenuation appears to be roughly constant 
in these bright galaxies, with slightly lower values of E(B- 
V) = 0.06 - 0.()8. We find that our median stellar masses 
woul d change by at most 0.14 dex (consistent with results 
from iPapovich et ^120011 that the choice of attenuation curve 
does not significantly affect stellar mass measurements), re¬ 
sulting in a minimal change to our major conclusion below 
on the stellar baryon fraction (although the slope of the stel¬ 
lar baryon fraction evolution assuming SMC dust is less than 
our fiducial scenario, and thus evolution is detected at reduced 
significance due to a lower value and larger uncertainty on the 
median mass of the z=7 sample in this scenario). We conclude 
that our assumption of a Calzetti et al. attenuation curve does 
not result in a strong bias in our results. 

The youngish character of these galaxies is surprising. Al¬ 
though age is notoriously difficult to measure robustly, the 
addition of the IRAC photometry does help, in particular at 
z = 6 , where the ensemble is constrained to have a typical 
age < 50 Myr ( 68 % C.L.). While evolved galaxies are not 
absent at these high redshifts (for example, z7_GNW_17001 
has log(M*/M 0 ) = 10.7 and an age of 400 Myr), they seem to 
be exceptions. This is in stark contrast with galaxies at z < 3 , 
where those residing at the bright-end of the luminosity func- 
ti on tend to be more e volved. 

iReddv et al.l (120061) published a stellar population analysis 
of galaxies at redshifts z 1-3.5. Using the results from their 
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models that assume a constant SFH, the median galaxy with 
4/uv < -21 (derived from the observed R-band magnitude 
and spectroscopic redshift) at z ^ 2 - 3 has log(M*/M 0 ) = 
10.2, an age of 260 Myr, E(B-V) = 0.16 and SFR = 70 
Mq yr“' (the results are similar when the sample is split into 
two redshift bins centered at z = 2 and z = 3). Thus, bright 
galaxies at the peak of cosmic star-formation activity have 
similar SFRs and dust attenuations as bright galaxies in the 
epoch of reionization, but the lower redshift galaxies are ~ 2 - 
2.5 X more massive and have rest-UV/optical light dominated 
by stars >5x older. Such a comparison with lower redshift 
makes sense if one is interested in the evolution of properties 
of similarly bright galaxies with redshift. It is also valid if we 
wish to explore the evolution of galaxies with Muv < Kv 
with redsh ift, as M,*,y at z 2.3 and 3 (-20.7 and -20.97, re¬ 
spectively ;|Redd^^[Steidel| 200 l) are similar to what we find 
at z > 4. 

A more interesting question is how UV-bright galaxies at 
z ^ 2-3 are related to UV-bright galaxies at z ^ 6 - 8 . Specif¬ 
ically, in view of the hierarchical galaxy assembly, are the 
former galaxies descendants of the latter? Are the latter 
progenitors of the former? Moreover, galaxy merging com¬ 
plicates direct number-counting-based matching across red- 
shifts. A few recent studies have tried comparing galax¬ 
ies at differ ent redshifts at the same cumulative number 
density (e.g.- Ivan Dokkum et alJ l2M^ iPanovich et al]l201 lb 

iLeia et alJl2013l) . iBehroozi et alJ (l2013ah showed that such a 

comparison is adequate for identifying th e low-z descendan t 
population of a high-z population (see also lJaacks et al.l2015[) . 
This is because the majority of massive high-z galaxies do 
not end up merging into substantially more massive systems. 
However, the converse is not true; reflecting hierarchal merg¬ 
ing, the cumulative number density of the high-z progenitor 
population of a low-z population has a comparatively higher 
cumulative number density. 

As shown in the following subsection, galaxies at z = 6 
and 7 have «(Muv < -21) = 3.6 x 10“^ Mpc“^ and 2.4 x 
10~^ Mpc~^, respective ly. Using the luminosity functions of 
iReddv & Steideil (I2009h atl.9<z<2.7 (the median redshift 
of the z ^ 2-3 comparison sample we quote here is z = 2.44), 
we And «(Muv < -21) = 2.24 x 10“"^ Mpc“^, unsurprisingly 
more abundant than our high-redshift sample. To match the 
cumulative number densities of our high-redshift samples, we 
need to select objects with Myv < -20.0 and < -19.8 for our 
z = 6 and 7 sample, respectively. The Muv < -21 galaxies at 
z 2-3 are thus plausible descendants of the high-z galaxies 
at these fainter magnitudes. 

Although above we only report stellar population results for 
brigh t galaxies, we performed SED fitting on the entire sam¬ 
ple of iFinkelstein et af] (120151) . finding that Muv = -20 galax¬ 
ies at z = 6 have mean stellar masses of log(M*/MQ) = 8.7 
and mean ages of 34 Myr. At z = 7, Muv = -19.8 galaxies 
have median stellar masses of log(M*/MQ) = 8.5 and median 
ages of 28 Myr. We conclude that, compared to either sim¬ 
ilarly bright or similarly abundant galaxies at lower-redshift, 
UV-bright galaxies at z > 6 are significantly younger and less 
massive, but have similar SFRs, as inferred from their similar 
U V luminosities, and exhibit relatively little evolution of dust 
attenuation at these luminosities. 

We acknowledge that ages are notoriously difficult to con¬ 
strain, as they are tied to the a ssumed star-formation histo¬ 
ries (e.g., iPapovich et alj|200l1) . In more realistic scenarios, 
the young stars that dominate the observed SED are mod¬ 


estly "outshining" the older generations, potentially biasing 
the measured age. However, the stellar masses (which we use 
for our primar y result in the fo llowing sections) are robust to 
age variations (iLee et al.ll201(lh . thus this potential bias does 
not affect our main results. 

3. HALO MASSES OF UV-BRIGHT GALAXIES 
3.1. Halo Masses 

Given the relatively high stellar masses and SERs seen in 
our sample of z > 6 galaxies, one may wonder if the processes 
governing gas cooling and the conversion of cold gas into stars 
differ from those at lower redshift. As a step toward answer¬ 
ing this question, we compare the stellar masses M* of these 
systems to the baryon masses in halos (Ub/Uni)Mh computed 
assuming a cosmic-m ean baryon fraction of Ub/Um = 0.1669 
(iKomatsu et al.ll 2 ()Tl]) . The clustering of these systems would 
have provided the most direct constraints on the halo mass, 
but the numbers and surface densities of these galaxies, partic¬ 
ularly at z > 7, are not y et sufficient to permit a robu st cluster¬ 
ing analysis (though see lBarone-Nugent et al.ll20T^ . Instead, 
we use abundance matching to estimate the halo masses for 
our bright galaxies. We refer the reader to previ ous works for 
a full discussion of abunda nce matching (e.g., iMoster et alJ 
l2010t IBehroozi et d1l2010h . Here, we review the procedure 
only briefly. 

Abundance matching assumes that the galaxy luminosity or 
stellar mass is a monotonic function of the halo mass. The 
most luminous galaxies are assumed to live in the most mas¬ 
sive halos. This is certainly a plausible assumption among the 
luminous galaxies we study here, but perhaps not among the 
more stochastic dwarf galaxies. Starting with a cosmological 
simulation (in this case, the pure ACDM simulation Bolshoi), 
one selects simulation snapshots close in redshift to the target 
redshift of an observational survey. In these snapshots, one 
selects a simulation volume equal to the volume of the sur¬ 
vey. One then identifies all virialized halos in the volume and 
rank-orders them by mass. After rank-ordering the observed 
galaxies by their luminosities, one places each galaxy in the 
simulated halo of the matching rank. This procedure provides 
a mapping of galaxy statistics onto host halo statistics. 

Here, we use Schechter parameterizatio n of the observed 
lumino sity functions at each redshift from lEinkelstein et alJ 
(120151) to estimate the host halo masses for the bright (M 1500 < 
-21) galaxies in our sample. Giyen the present stellar mass 
function uncertainties, luminosity-based matching is cur¬ 
rently more robust than the stellar mass-based matching. In 
the future, the matching should be performed directly with 
the stellar mass as it should be more correl ated than the UV 
luminosity with the host halo mass (e.g., ILee et al.l IlM 
iGerke et aUlIoil) . 

Eor our analysis, we use the results from the Bolshoi cos¬ 
mological simulations (iKlypin et alJl201 1 I) . which has 2048^ 
particles in a 250 (Mpc//z)^ box. This translates to a halo 
mass resolution of log(Mh/MQ) = 10. As we will find that our 
galaxies haye halo masses ^ 10'*’ Mq, the resolution of Bol¬ 
shoi is more than sufficient for the present analysis . We us e 
the halo mass functions deri yed in IBehroozi et al.l (I20131J) . 
which are a modification of the lTinker et al.l(l2008h mass func¬ 
tions, include subhalos, and are accurate at ye ry high red- 
shifts (see Appendix G of iBehroozi et al.ir2013bl) . We deriyed 
halo mass functions at our redshifts of interest by yolume- 
ayeraging the Bolshoi snapshot mass functions oyer the same 
redshift ranges as those defining our galaxy samples. 
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Fig. 3.— Top left: The cumulative luminosity function at z = 4, 5, 6, and 7. Bottom Left: Cumulative halo mass functions at z = 4, 5, 6, and 7, derived by 
volume-averaging the Bolshoi snapshot mass functions over the same redshift ranges as those defining our galaxy samples. The arrows show our results from 
abundance matching at z = 7, where galaxies with M 1500 < “21, which have «(Mqv < —21) = 2.5 x 10“^ Mpc^, have halo masses of log(Mh/MQ) = 11.35. Right: 
Relation between observed UV absolute magnitude and abundance-matching-derived halo mass at our redshifts of interest. The arrows denote the halo masses at 
our magnitude of interest of Muv = -21. 


The top-left panel of Figure [3 shows the cumulative lumi¬ 
nosity functions at z = 4-7 and the bottom-left panel shows the 
cumulative halo mass functions at the same redshifts. To infer 
the host halo masses for the bright galaxies with Misoo < -21, 
we first find the cumulative number density of galaxies at 
that magnitude. For example, for z = 7, the number density 
is indicated with the horizontal arrow in the two left panels 
of Figure [3 The host halo masses can then be inferred by 
finding the halo mass above which the cumulative halo num¬ 
ber density equals the cumulative galaxy density. We find 
log(Mh/MQ)= 11.93, 11.68, 11.57, and 11.35 at z = 4, 5, 6, 
and 7, respectively. To estimate the uncertainties in the halo 
masses, we ran 10^ Monte Carlo simulations, in each drawing 
a luminosity function randomly from the MCMC sample gen¬ 
erated during the luminosity function estimation, which ac¬ 
count for both Poisson noise and uncertainti es in the luminos¬ 
ity fu nction completeness simulations (see iFinkelstein et al.l 
120151) . We find that the uncertainties in the halo mass are low, 
~0.01-0.03 dex, reflecting relatively low uncertainties in the 
cumulative luminosity functions. 

3.1.1. Cosmic Variance 

To estimate how our measured abundances of bright galax¬ 
ies are affected by cosmic variance, we use d a suite of semi- 
analyt ic models. These models, based on ISomerville et alJ 
(120081) . were provided to the CANDELS team, and include 
a set of mock catalogs, one for each of the five CANDELS 
fields (with each individual mock catalog covering a volume 
somewhat larger than the observed volume). We tuned the 
dust attenuation in these SAMs to match the obse rved UV lu¬ 
minosity functions from lEinkelstein et alJ (l2015h . In each of 
our redshift bins, we used the SAMS to extract 64 indepen¬ 
dent volume elements comparable to one GOODS-sized field 
(16' X 10') and 3192 independent volume elements compa¬ 
rable to a single WEC3 pointing (2.1' x 2.1'). In each red¬ 
shift bin, the fractional uncertainty due to the combination of 
Poisson fluctuations and cosmic variance was derived as the 
standard deviation (computed over all of the realizations of a 
given field size) of the number of galaxies in a luminosity bin 


of AMuy = 0.5 mag centered at Myy = -21, divided by the 
mean number of galaxies in the bin. 

The total fractional uncertainty due to Poisson fluctuations 
and cosmic variance was then derived by combining the vari¬ 
ances from two GOODS-sized fields and five single WPC3 
pointing-sized fields, as th is was the area used in th e luminos¬ 
ity function calculation bv lEinkelstein et aH (12015h . who used 
GOODS-S, GOODS-N, and five individual fields (the HUDE, 
the two HUDE parallels, and the first two first-year Erontier 
Eields parallel fields). We find that the fractional uncertainties 
due to cosmic variance are 0.132, 0.159, 0.212, and 0.327 at 
z = 4, 5, 6, and 7, respectively. We include these uncertainties 
in the Monte Carlo simulations discussed in the above para¬ 
graph, and find that the uncertainty in the derived halo mass 
increases by a factor of ^ 2 to 0.03 dex at z = 4 and 0.06 dex at 
z = 7 with the inclusion of the uncertainty due to cosmic vari¬ 
ance. Our derived halo masses, and these total uncertainties, 
are listed in Table 2. 

3.2. Evolution of the Stellar-to-Halo Mass Ratio 

Comparing the halo masses estimated in the previous sub¬ 
section to the median stellar masses estimated in §2.3, we can 
calculate the ratio of the median stellar mass to halo mass 
(SMHM). We find M*/Mh = 0.009,0.013, 0.016, and 0.020 at 
z = 4, 5, 6, and 7, respectively. Thus, at a constant UV lumi¬ 
nosity, the stellar-to-halo mass ratio increases with increasing 
redshift. These results are listed in Table 2. 

The cosmic baryon mass fraction is much higher than the 
MjjM-a ratios we derive, at Ub/Hnj = 0.1669 (iKomatsu et alJ 
120111) . In EigurelH we show the evolution of the SMHM ra¬ 
tio in the units of Hb/Om as a function of redshift. We re¬ 
fer to this quantity as the stellar baryon fraction (SEE), as it 
measures the amount of baryons converted into stars com¬ 
pared to the cosmic allotment of baryons in the halo. We 
find that the SEE evolves from 0.117 ±0.043 at z = 7 to 
0.051 ± 0.006 at z = 4. This factor of ^ 3 evolution is sig¬ 
nificant, as fitting a linear function for SEE(z) yields the slope 
c/SEE/c/z = 0.0239 ± 0.0074, with the trend detected at 3.2a 
significance. 
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Fig. 4.— The stellar baryon fraction (SBF) in bright (M\]\ = —21) galaxies from z = 4 to 7. We define the SBF as the stellar to halo mass ratio in units of the 
cosmic baryon mass fraction Ob/fim. We find that the SBF increases with increasing redshift, which may be responsible for the apparent lack of evolution in the 
characteristic magnitude observed over this redshift range. 


The inferred significance of the trend depends on our as¬ 
sumed uncertainties in the stellar baryon fractions. We as¬ 
sumed an uncertainty in the halo mass and median stellar mass 
as reported in Table 2. The quoted stellar mass uncertainties, 
which decrease from 0.14 dex at z = 7 to 0.04 dex at z = 4, are 
well below the uncertainty for an individual object, as here we 
are interested in the accuracy with which we can constrain the 
median stellar mass at different redshifts. However, these val¬ 
ues are consistent with the typical u ncertainty in the m edian 
stellar mass at Myv = -21 derived bv ISong et al.l (12015h . who 
fit the mass-to-light ratio over a wide dynamic range in UV 
luminosity, finding that it decreases from 0.1 dex at z = 7 to 
0.02 dex at z = 4. If one used the median of individual galaxy 
mass uncertainties in each of our samples, which decreases 
from 0.19 dex at z = 7 to 0.11 dex at z = 4, the trend of in¬ 
creasing SBF with redshift is still apparent (albeit reduced in 
significance to ~ 1.7 ct). 

Finally, we recall that our fiducial sample was selected to 
include galaxies with Muv < -21. Given the shape of the 
luminosity function, the median magnitude of this sample is 
very close to -21 (ranging from -21.2 to -21.3). To check 
if our sample selection biases the median masses, we exam¬ 
ined how the evolution of the stellar baryon fraction changes 
if we use the median stellar mass of galaxies with luminosi¬ 
ties Muv = -21 ±0.25. We find a comparable number of 
galaxies as in our fiducial sample, specifically 118, 76, 20, 
and 17 galaxies at z = 4, 5, 6 , and 7, respectively. The me¬ 
dian stellar mass is slightly lower than in our fiducial sam¬ 
ple, log(M*/M 0 ) = 9.68, 9.66, 9.50 and 9.52 and z = 4, 5, 
6 and 7, respectively. The amplitude of the redshift deriva¬ 
tive of the stellar baryon fraction is thus lower, c/SBF/c/z = 
0.0123 ± 0.0041, yet the evolution is still significant at the 
3.1cr level. We conclude that there is significant evolution in 
the stellar baryon fraction in that it decreases with decreasing 
redshift, and that this evolution is stable against several def¬ 


initions of both the median stellar mass and the stellar mass 
uncertainty. 

3.2.1. Comparison to Previous Results 

iBehroozi et al.l (120131j) recently studied the evolution of 
the SMHM relation by modeling all available observational 
constraints, including luminosity functions, stellar mass func¬ 
tions, and SFRs, and exploring the galaxy evolution parame¬ 
ter space with an MCMC search. They found that the SMHM 
curve peaks at a roughly constant halo mass of log(Mh/MQ) = 
11.7 at z = 0^. Specifically, at z = 4, they found a peak 
SBF at a halo mass of log(Mh/M 0 ) Ri 12.0, consistent with 
the typical halo mass we derive for our bright z = 4 galax¬ 
ies, log(Mh/M 0 ) = 11.9. Then they went on to find that the 
peak of the SMHM relation at z > 5 occurs at a halo mass 
steadily decreasing with redshift, log(Mh/M 0 ) = 11.9, 11.6, 
and 11.4 at z = 5, 6 , and 7, respectively. This is very similar 
to what we find for the typical halo masses of M^y galaxies, 
log(Mh/M 0 ) = 11.7, 11.6, and 11.4 at z = 5, 6 , and 7, respec¬ 
tively. The peak SMHM ratio measured by IBehroozi et alJ 
(120131j) is somewhat higher, ^^3.7% at z = 7 (converting from 
a Chabrier to our Salpeter IMF) than the 2.0% we find here. 
However the two values are indistinguishable given the sig¬ 
nificant uncertainties in the high-redshift ob servables, par¬ 
ticularly the stellar mass functions used by IBehroozi et alJ 
(I2013bl) . 

Although here we are specifically concerned with the halo 
masses of bright galaxies, in the right-hand panel of Figure [3] 
we also provide the abundance-matching-derived halo masses 
at all observed magnitudes in each of our redshift bins. Our 
derived halo mass of log(Mh/M 0 ) = 11.35 at z = 7 is consis¬ 
tent with a recent clustering-based mea sure of log(Mh/M 0 ) k, 
11.2 from iBarone-Nugent etalJ (120141) . However, their halo 
mass estimate was for galaxies with Muv < -19.4, a sam¬ 
ple that has a lower average luminosity than our sample, and 
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TABLE 2 

Dark Matter Halo Properties of Bright Galaxies 


Redshift 

log n (Muv < -21) 

(Mpc“^) 

logMh 

(Mq) 

Median 

M*/Mh 

Stellar Baryon 
Fraction 

z = 4 

-3.86 

11.93+^“ 

0.009±0.001 

0.051 ±0.006 

z = 5 

-4.01 

11.68|i 

0.013±0.002 

0.079±0.013 

z = 6 

-4.45 

11.57+«-“ 

0.016±0.003 

0.097±0.019 

z = 7 

-4.62 

11 35® 
ii.jo_0 06 

0.020±0.007 

0.117±0.043 


Note. — The uncertainties in the halo mass are derived via Monte 
Carlo simulations and include the uncertainty in the number density of 
Muv < -21 galaxies, which reflects our fiducial luminosity function un¬ 
certainties as well as cosmic variance. The uncertainties in M*/Mh and 
the stellar baryon fraction assume an uncertainty in the median stellar 
mass from Table 1. 


thus a lower halo mass is expected, as shown in the right-hand 
panel of Figure 3. Our z = 4 re sults are consiste nt with earlier 
clustering-based estimates by iLee et al.1 (|200^, while at z = 
5 and z = 6, our derived halo masses are somewhat hi gher 
than the clustering-ba sed estimates from iLee et al.l (I2006^ and 
lOverzier et al.l (l2006h . respectively, likely due to the fainter 
luminosities considered in those works (Muv < -20 -19.5, 
respectively). Although there are minor differences, it is en¬ 
couraging that two independent methods, abundance match¬ 
ing and clustering, tentatively agree on the halo mass esti¬ 
mates at such high redshifts, though certainly the clustering- 
based can be made more robust at z > 6. 


3.2.2. UVLuminosity Scatter 

The halo mass estimates in § 13.11 did not explicitly model 
the effect of the scatter in UV luminosity at a fixed halo 
mass. The estimates compared the number density in galaxies 
with Muv < -21 to the number density in halos with masses 
> Mh to derive the characteristic halo masses Mh. The halo 
mass-luminosity relation has an intrinsic scatter that is usually 
treated as a Gaussian in the magnitude Muv centered on some 
median (or mean) halo-mass-dependent magnitude. Because 
of this scatter, some halos with masses > Mh can can host 
galaxies with atypically low luminosities, Muv > -2 1, p lac- 
ing them outside of our luminosity cut. However, in § 13.11 the 
halos hosting these faint interlopers were being counted in the 
cumulative halo mass function and this could have biased the 
derived Mh. If we had known which halos hosted faint inter¬ 
lopers (which we do not), we would have excluded them from 
the halo counting to derive a more accurate Mh. Excluding a 
fraction of halos at every halo mass, the resulting characteris¬ 
tic halo mass Mh would have b een lower than the scatter-blind 
estimate (iBehroozi et al.ll20T^ . 

Here, we attempt to quantify this bias. We estimate the 
amplitude of the luminosity scatter from the relati on between 
the UV luminosity and the stellar mass derived bv ISong et al.l 
(HoBl). This is warranted because the stellar mass is ex¬ 
pected to be more correlated with the halo mass than with 
UV luminosity. We find that at the mass of our sample 
log(M*/MQ) ^ 9.8, the scatter in the UV absolute magnitude 
is ^0.3 dex at all redshifts we consider. 

To assess the impact of luminosity scatter on our abundance 
matching-determined halo mas ses, we carried out t he itera¬ 
tive deconvolution described in iReddick et al.l (120131) to get a 
handle on the true, underlying bivariate distribution of halo 
masses and observed UV luminosities that exhibits an intrin¬ 
sic luminosity scatter at fixed halo mass. We performed abun¬ 


dance matching on the bivariate distribution properly apply¬ 
ing the Muv < -21 cut to find halo masses log(Mh/MQ) = 
11.65, 11.44, 11.36, and 11.13 at z = 4, 5, 6, and 7, respec¬ 
tively. These values are ~ 0.2 - 0.25 dex lower than the val- 
ues obtained by scatter-blind abundance matching in § ITT] 
Therefore we expect that the true host halo masses of our 
Muv < -21 galaxies are slightly lower than those reported in 
Table 2. If the luminosity scatter had been 0.5 dex, our halo 
masses would be overestimated by 0.7-0.9 dex. 

Since halo masses at all redshifts are corrected by a similar 
logarithmic increment, the effect of luminosity scatter does 
not affect our detection of an evolving stellar baryon fraction. 
But for higher luminosity scatter which could characterize a 
different regime in which galaxy formation is more stochastic, 
this effect could be more significant. 

3.2.3. Systematic Uncertainties in Stellar Mass 

Our derivation of the stellar masses required a variety of 
assumptions that could have potentially biased our results. 
Could any of the biases explain away our conclusion that UV- 
bright galaxies from z = 4 to 7 have similar stellar masses but 
are found in progressively lower mass halos toward higher 
redshifts? One strong assumption we made is that of the 
Salpeter IMF. If the high-mass slope of the IMF evolved 
with redshift, this certainly could have biased the conclusion. 
Other assumptions involve the parameterization of SFHs and 
the treatment of metallicities and dust attenuations. As dis¬ 
cussed earlier, the UV-continuum slopes in UV-bright galax¬ 
ies appear to be roughly constant across this redshift range, 
and so dust abundances will only affect our results if the typ¬ 
ical dust-law changes as a function of redshift. Direct mea¬ 
surement of the metallicities is be yond our current capabi lities 
at these high redshifts (though see lFinkelstein et al.l2013[) . but 
a plausible metallicity variation produces only a minor change 
in colors. The SFH parametrization is potentially more trou¬ 
bling, though recent results indicate that at least on aver- 
age, galaxy-scale S FRs are smooth functions of time (e.g., 
ISalmon et al.ll2015h . 

A basic test of any bias in the mass measurements is to com¬ 
pare the shapes of the SEDs of our galaxies. Eigure|5]shows a 
median flux stack of galaxies in each of our redshift bins ver¬ 
sus the rest-frame wavelength, scaled vertically to a common 
redshift z = 6. The shapes of the SEDs are remarkably simi¬ 
lar, especially in the rest-frame UV, where they appear iden¬ 
tical. Modest differences are visible in the rest-frame optical, 
likely due to the lower signal-to-noise of Sp/fzer/IRAC data, 
as well as the strong nebular [O III] and Ha lines which red¬ 
shift through the bandpasses. Given the highly similar SED 
shapes, it is unlikely that an unaccounted for systematic ef¬ 
fects strongly bias our results. Rather, we appear to be study¬ 
ing a very similar type of galaxy at each redshift; this type 
of galaxy lives in lower mass halos at higher redshift. This 
conclusion is confirmed by stellar population model fits to the 
stacks, which yield stellar masses consistent within ^0.3 dex 
of the median stellar masses in Table 1 (log[M*/MQ] = 9.8, 
10.0, 9.9 and 9.9 at z = 4, 5, 6, and 7, respectively). 

3.2.4. Dusty Star-Forming Galaxies 

Our sample is selected on the basis of UV luminosity and 
therefore it is prudent to examine what effect that may have 
on our results. In particular, a rest-frame UV selection may 
miss extremely dusty galaxies which have their UV light at¬ 
tenuated below our detection sensitivity. Such systems, re- 
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Fig. 5.— The median flux-stacked SEDs of our samples of bright galaxies 
at each redshift. The spectral shapes of these SEDs are remarkably similar, 
with minor differences appearing at longer wavelengths due to the presence 
of nebulai' emission lines, as well as the generally lower signal-to-noise of 
the long-wavelength data. The near identical nature of these SEDs confirms 
our conclusion that bright galaxies from z = 4 to 7 are physically very similar, 
and that these log(M* /Mq) = 9.6-9.9 galaxies on average inhabit lower mass 
halos at higher redshifts. 

ferred to as sub-millimeter galaxies (SMGs; after their se¬ 
lection wavelength), have been well studied at moderate red- 
shifts. Thes£_galaxie^re typically very massive (log M/Mq 
> 11: ICasev et al.ll20ri) . If we were missing a large popula¬ 
tion of these galaxies, it would bias our derived halo masses to 
be too high (as we would be placing our observed galaxies in 
the most massive halos in our volume, which would truly be 
occupied by these dusty galaxies). If the abundance of SMGs 
evolves with redshift at z > 4, then not accounting for these 
systems would bias our halo masses differently at different 
redshifts, biasing our inferred evolution of the stellar baryon 
fraction. 

To explore the potential bias introduced by our UV-only se¬ 
lection, we compare the space density of SMGs at z = 4 to our 
UV luminosity function. The space density of SMGs at high 
redshift is very uncertain for a variety of reasons, including 
the relatively shallow depths of sub-millimeter surveys, and 
the difficulty of obtaining redsh i fts for such systems. We start 
at z = 2, where IChapman et al.l (l2005h found a space density 
for SMGs of 5 X 10 ^ Mpc A more complet e survey for 
SMGs at multiple wavelengths was performed bv ICasev et al.l 
(120131). who found a SMG SFR density ^2x higher (see Fig. 
23 of ICasev et alJ 120141) . implying a space density of z = 2 
SMGs of 10“"* Mpc To estimate the evolution in this 
quantity to z = 4, we use the redshift evolution shown by 
ICasev et’aD (l2014h . which shows that the SFR density from 
SMGs is an order of magnitude lower at z ~ 4 than at z 2. 
Therefore we adopt 10“^ Mpc as a fiducial space density 
forz = 4SMGs. 

We simulated the absence of SMGs in our sample by adding 
this abundance to our observed cumulative luminosity func¬ 
tions, and repeating the abundance matching analysis. The 
absence of such a population of SMGs results in a bias of the 
z = 4 halo mass for our galaxies of interest of only 0.02 dex 
(11.91 versus our fiducial result of 11.93). This is within the 
68% confidence range on our fiducial halo mass, therefore not 
significant. If the abundance of SMGs stays constant to z = 
7, the declining UV luminosity function results in a slightly 
larger overestimate of the halo mass of 0.06 dex at z = 7. Were 


this the case, the observed evolution in the stellar baryon frac¬ 
tion would be even stronger than we observe. However, such 
a high abundance of dusty galaxies at z = 7 is highly unlikely. 
We therefore conclude that our selection does not affect our 
main conclusions in this study. 

4. DISCUSSION 

To understand the physical effects responsible for our ob¬ 
served trend of an increasing SBF with increasing redshift, 
here we consider a variety of possible mechanisms, with the 
caveat that our observations cannot uniquely distinguish be¬ 
tween these scenarios. First, galactic gas at higher redshifts 

has higher surface densities Sg oc /gasA/j, (1 -Fz) , where /gas 
is the gas fraction in the cooled, virialized phase. The typi¬ 
cal free-fall time to which the gas-to-stars conversion rate is 
proportional (albeit with a small coefficient—dimensionless 
SFR—see, e.g., iKrumholz et al.l |2012^ is fpp oc (1-Fz)“^/^. 
However, masses of the most massive progenitors (MMP) of 
our galaxies’ host halos are Mh,MMp(zO = Mh( z)e~°^^ with 
a K, 1 where z' > z is the p rogenitor redshift (iNeistein et alJ 
l2006t iFakhouri et alJ 120101) . Therefore the mass-doubling 
growth time is fgrow(z) ~ [(1 +z)H{z)d\nM'a_yaAp/dzT'^ oc (1 -F 
z)-5/2 jjjeans that the ratio of the free-fall time to the 

growth time increases with increasing redshift as fpp/fgrow oc 
1 -F z, suggesting that if the minor progenitors are inefficient 
star formers so that they do not contribute substantial stel¬ 
lar mass to the main branch, it should be progressively more 
difficult at high red shifts to convert the g as acquired through 
growth into stars (iBouche et al.l 120101 : IKrumholz & Dek^ 
1201 2t iDekel et^ 120131) . Second, the gas-phase metallicity 
and associated dust abundance appear to decrease with in¬ 
creasing redshift and decrea sing mass ('e.g.. lFinkelstein et alJ 
l2012btlBouwens et al.l2014l) . Because metals and dust are the 
principal gas-cooling and UV-shielding agents, the observed 
trend of decreasing dust with increasing redshift could have a 
dramatic effect on the abundance of the cold (T < lOOOK) 
gas in which star formation seems to exclusiv ely happen 
(IKrumholz & Dek^ 120121 IKrumholz et al.ll^l2^ . The drop 
in cold gas abundance toward higher redshifts would imply 
lower SFRs. 

All this suggests that from the supply-versus-consumption 
side alone, we expect an opposite trend from the one mea¬ 
sured, a trend in which the stellar baryon fraction decreases 
with increasing redshift. How should we then interpret the 
stellar baryon fractions that increase with increasing redshift? 
We can only speculate. Our observations imply that cold gas 
is less readily available for star-formation at lower redshifts. 
How might this come about? One scenario involves the tran¬ 
sition of gas from the warm to cold phase, which occurs at 
a density ~ 0.1-lcm“^ (e.g.. IWolfire et al.]ll995L l2003h . As 
the gas density decreases towards lower redshift, a progres¬ 
sively smaller mass fraction of the neutral gas phase is cold. 
Additionally, the nature of gas collapse in the disk is criti¬ 
cal, as the galaxy-wide star formation rate may be particu¬ 
larly sensitive to whether the central gaseous str ucture (typi¬ 
cally a clumpy disk) is viol ently self-gravitating (iDekel et alJ 
120091: ICeverino et al.l 120101) . Star formation is most efficient 
when large gas clumps become self-gravitating, and large 
self-gravitating clumps of gas have b een observed to fo rm 
stars with high efficiency at z ~ 3 (e.g.. lGenzel et alJl201 ll) . 

Another possibility is that we are witnessing the effect 
of the growth of the circumgalactic medium (CGM) around 
galaxies. Observations with the Cosmic Origins Spectrograph 
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on HSThwe, recently revealed circumgalactic ionized gas ha¬ 
los at low redshift con taining significant baryon and synthe¬ 
sized metal fractions (IWerk et al.l 120141) . The CGM likely 
grew over time, as shock waves driven by supernova blast 
waves and superbubbles raised an increasing cumulative gas 
fraction to high temperatures. This reprocessed gas would be 
deposited in a warm ionized CGM and, at least for a period of 
time, off limits to conversion into the cold phase. However, 
the physics of the CGM is poorly understood, and in particu¬ 
lar it is not clear at what rate the CGM gas recondenses back 
into the warm neutral phase in the disk. 

Finally, the feedback which builds the CGM may be less ef¬ 
ficient at high redshift, further increasing the amount of cold 
gas available for star formation. One may expect that the ef¬ 
fects of feedback are less significant at high redshift where 
halos are denser and at fixed mass have higher circular ve¬ 
locities. Both make it more difficult for SNe to power galac¬ 
tic outflows that eject material into the CGM or outs i de the 
halo. Using the virial estimate from iBrvan & NormanI (119981) 
we calculated the circular velocities of the host halos in our 
galaxy sample to be 242, 219, 219 and 197 km s“* at z = 
4, 5, 6 and 7, respectively. Thus, the circular velocities in 
fact decrease towards higher redshift, which would if any¬ 
thing make it easier to eject material at higher redshift. How¬ 
ever, the properties of the feedback mechanisms may evolve 
as well. First, as galaxies evolve and the typical gas density 
decreases, the supernova remnant thermalization efficiency, 
which is an increasing function of the cooling time in HII re¬ 
gions shock-heated by supernova blastwaves, increases, i.e., a 
progressively smaller fraction of the ini tial 10 ^* erg per su - 
pernova is quickly irradiated (see, e.g.. lCreasev et al.l 120131) . 
This implies that at lower redshifts, a larger fraction of super¬ 
nova energy may be available to drive galaxy-wide outflows, 
thus making the feedback from star formation more effective. 
It is also possible that due to the very short growth times at 
high redshifts, so much gas is piling on that the outflows are 
somehow bottled in. Furthermore, the higher dust content at 
lower redshift can lead to stronger momentum-driven radia¬ 
tive stellar feedback, regulating further star formatio n (see, 
e.g.. lMurrav et al.ll2010l : Andrews & Thomr)sonll20lll) . 

One additional potential physical mechanism which may 
evolve is the ability of feedback from active galactic nuclei 
(AGNs) to suppress star formation, which is commonly im¬ 
plemented in theoretical models to avoid an overabundance of 
brigh t/massive galaxies (see discussion in ISomerville & Davd 
|203)- This type of feedback requires an accreting super- 
massive black hole, and although there are some examples 
of bright AGNs at very high redshift, the AGN/quasar lu- 
minosity function ap p ears to decrease rap idly at z > 3 (e.g., 
iRichards et al.ll200^ . iBowler et al.l (12015h have recently ob¬ 
served that the bright end of the galaxy UV luminosity func¬ 
tion was steeper than the halo mass function at z = 6 , but not 
at z = 7, and hypothesized that such an observation could 
be explained if feedback in bright/massive galaxies due to 
AGN first “turned on” at z < 6 . However, the details of how 
AGN couple with galaxies and their surroundings, particu¬ 
larly at these epochs, are highly uncertain, so it remains un¬ 
clear whether black hole accretion has significantly affected 
the growth of the galaxies we consider here. 

While the scenarios we have discussed are clearly specula¬ 
tive, our observations imply that the latter effects, primarily 
a reduced efficiency of feedback at higher redshift due to a 
variety of redshift-dependent effects, control the evolution of 


the stellar baryon fraction. We conclude that the true cause of 
how a larger fraction of the baryons turns into stars at higher 
redshifts is most certainly determined by a delicate competi¬ 
tion of factors. 

Regardless of the underlying cause, the consequences of the 
increased availability of cold gas are intriguing. Our results 
at z = 6 and 7 show that ~10%-12% of the cosmic comple¬ 
ment of baryons in these galaxies has been converted into 
stars. The remaining baryons must exist in the gas phase. 
If they are in the warm or cold neutral phase or the molec¬ 
ular phase, then the gas fract ion in these phases is much 
higher than at lower redshifts (iMagdis et al.l 1201^ . A high 
gas mass fraction at very high redshifts is not unexpected and 
may soon be confirmed the Atacama Large Millimeter Array 
(ALMA; e.g., by measuring the dynamical mass via spectrally 
resolved FIR emission lines). There are indirect empirical 
hints tha t neutral and molecula r gas fractions increase with 
redshift. iPapovich et alJ (1201 Ih studied the evolution of ob¬ 
served galaxy SFRs, stellar masses, and sizes, and concluded 
that the gas-to-stellar mass fraction must rise with redshift, 
reaching Mgas/M* = 3.9 at z = 7 (at a higher cumulative num¬ 
ber density of 2 x 10“^ Mpc“^). This could intriguingly play 


some role in the decreasing visibilip 

of Lya at z > 6 (e.g.. 

Fontana et alJ 120101: 

Pentericci et alJ 

20l1tlTi1vi et alJl20l4 

Pentericci et alJl2014 

). 


Finally, we examine the likely descendants of these bright 
z = 4-7 galaxies. Specifically, being among the most rapidly 
star forming galaxies at their redshifts, could some of them 
end up evolving into extreme systems such as sub-millimeter 
galaxies (SMGs) b y lower redshifts? Us ing the halo mass 
evolution tool from iBehroozi et aTI (l2013al) . we calculated the 
68 % confidence range of the descendant halo masses of our 
sample of bright high-redshift galaxies. We And that by z = 2, 
galaxies we observe at z = 4, 5, 6 , and 7 will exist in halos with 
log(Mh/Mo) = 12.3-12.8, 12.3-12.8, 12.4-13, and 12.4- 
13.1, respectively. SMGs ar e thought to be hosted by halos 
with l og(Mti w 13 (e.g.. lHickox et al.ll2012tl^asev et alJ 
l2014h . thus the majority of M^y galaxies at z = 4-5 will evolve 
into galaxies at z = 2 with halos slightly less massive than 
those of the typical SMGs. However, the SMG host halo mass 
begins to be consistent with the expected descendants of the 
z = 6-7 galaxies, thus some subset of very UV-bright galaxies 
at z > 6 may, in principle, turn into lower-redshift SMGs. A 
basic test of this is to see whether, if we assume these galaxies 
keep their current SFRs, they can grow large enough to match 
the stellar mass of a typical SMG by z = 2. If we assume a SFR 
of 50 Mq yr“* (Table 1), we find that the z = 2 descendants of 
Myy galaxies at z = 6 and 7 will have log(M*/M 0 ) = 11.1. 
This is in the range of SMG stellar mas ses found in the lit¬ 
erature (see review in ICasev et^l2014f) . Our assumed con¬ 
stant SFR of ^ 50Mq yr“*, which implies stellar masses con¬ 
sistent with those observed in SMGs, is approximately the 
SFR disk galaxies should have prior to coalescence to pro- 
duce SMGs in the merg er-driven scenario for SMG formation 
(iNaravanan et al.l[201(il) . 


5. CONCLUSIONS 

Recent observations have shown that the characteristic UV 
luminosity M^y does not significantly evolve from z = 4 to 7, 
which is unexpected given the general decline in the cosmic 
SFR density towards higher redshift over that time. To investi¬ 
gate the physical effects behind this observed non-evolution in 
Myy, we have inspected the stellar populations in M^y galax- 
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ies at z = 4 to 7. We have found the following results: 

• Galaxies with Myv < - 21 appear to have very similar 
physical properties at z = 4, 5, 6 , and 7, including stellar 
mass, dust attenuation and SFR. 

• Using abundance matching to infer the likely hosting 
halo masses, we found that Muv < -21 galaxies, which 
we have measured to have log(M*/M 0 ) =9.6-9.9, live 
in progressively smaller halos towards higher redshift, 
with log(Mh/MQ) = 11.93 at z = 4, and log(Mh/M 0 ) = 
11.35 atz=7. 

• The stellar baryon fraction, defined as the fraction of 
baryons in stars in units of the cosmic mean Ub/Um 
rises from 0.051 ±0.006 at z = 4 to 0.117 ±0.043 at 
z = 7. This trend is significant at the 3.2 (t level. 

• The observed trend of an increasing SBF with increas¬ 
ing redshift does not agree with simple expectations of 
how galaxies grow. Therefore, our observations im¬ 
ply a change in the physical properties governing star- 
formation at z > 4, such as, for example, a reduced effi¬ 
ciency of stellar and supernova feedback toward higher 
redshifts. 

Future studies can improve upon our results by probing a 
larger volume to increase the sample of bright galaxies, al¬ 
lowing us to establish the evolution of the SBF at greater 
significance, as well as extending this analysis to z = 8 and 
9. Additionally, a more robust determination of the ratio be¬ 
tween stellar mass and halo mass, and thus the SBF, can be 
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APPENDIX 

In the above text, we discuss the median stellar population properties of our sample of galaxies. Galaxies of course are an 
inhomogeneous population, thus in this appendix we provide tables of the individual properties for each of the bright galaxies 
used in our analysis (those with IRAC measurements, and a best-fitting SED model with < 20; see §2.3). Each column in 
these tables lists the best-fitting value along with the 68% confidence range in parentheses. 

TABLE Al 

Stellar Populations of Bright Galaxies at z = 4 


ID 

Right Ascension 
(J2000) 

Declination 

(J2000) 

Redshift 

logM* 

(Mq) 

Age 

(Myr) 

E(B-V) 

SFR 

(M 0 yr-') 

z4 GSD 34736 

53.096840 

-27.866074 

3.51 (3.44^3.66) 

9.97 (10.13-10.49) 

10(19-10!) 

0.34 (0.18-0.34) 

451 (109^34) 

z4 GSD 30292 

53.086891 

-27.844139 

3.51 (3.46-3.75) 

10.49(10.31-10.51) 

57 (49-90) 

0.24(0.16-0.32) 

132(64-322) 

z4 GNW 9013 

189.085114 

62.160465 

3.54 (3.49-3.60) 

9.97 (9.93-10.15) 

1015(202-1015) 

0 . 10 ( 0 . 02 - 0 . 10 ) 

54 (26-53) 

z4 ERS 4095 

53.143330 

-27.690090 

3.55 (3.48-3.67) 

10.27 (10.17-10.55) 

49 (30-101) 

0.16(0.14-0.28) 

88(73-299) 

z4 ERS 9969 

53.120926 

-27.709446 

3.57 (3.52-3.64) 

9.68 (9.43-9.68) 

57 (30-286) 

0.04(0.00-0.10) 

20(14-38) 

z4 GNW 26176 

189.483200 

62.284786 

3.63 (3.57-3.71) 

9.99 (9.96-10.18) 

90 (40-202) 

0.12(0.02-0.14) 

65 (27-85) 

z4 GNW 1986 

189.159058 

62.115471 

3.64 (3.58-3.70) 

9.98 (9.96-10.19) 

101(49-1015) 

0.08 (0.06-0.10) 

55 (44-77) 

z4 ERS 16929 

53.087231 

-27.729538 

3.64 (3.58-3.72) 

9.42 (9.21-9.85) 

19(10-40) 

0.14(0.14-0.20) 

59 (55-109) 

z4 GSD 29028 

53.087368 

-27.839535 

3.64 (3.52-3.75) 

10.17 (10.12-10.22) 

80(80-101) 

0.10(0.06-0.10) 

38 (25^0) 

z4 ERS 20075 

53.020580 

-27.742151 

3.65 (3.59-3.71) 

10.55 (10.54-10.66) 

202(202-570) 

0.06 (0.02-0.08) 

75(53-104) 

z4 GND 30689 

189.339355 

62.216450 

3.66 (3.59-3.75) 

9.45 (9.40-9.70) 

39 (30-71) 

0.04(0.00-0.10) 

20(14-36) 

z4 GSW 4356 

53.109478 

-27.879360 

3.67 (3.61-3.77) 

9.72 (9.66-9.90) 

57 (49-90) 

0.00 (0.00-0.08) 

19 (19^3) 

z4 GSD 15786 

53.071735 

-27.798437 

3.67 (3.60-3.75) 

9.61 (9.67-9.88) 

49(57-1015) 

0.04(0.00-0.10) 

19(14-33) 

z4 GSD 535 

53.198959 

-27.737940 

3.69 (3.57-3.81) 

9.59 (9.53-9.79) 

49(40-57) 

0.00 (0.00-0.08) 

18(17-39) 

z4 ERS 3396 

53.117710 

-27.686771 

3.69 (3.64-3.75) 

9.90 (9.53-9.90) 

49 (30-49) 

0.10(0.08-0.16) 

48 (40-74) 

z4 GSD 21002 

53.121414 

-27.814621 

3.69 (3.63-3.76) 

9.80 (9.57-9.85) 

30(10-30) 

0.18(0.16-0.24) 

93 (70-152) 

z4 GSD 35257 

53.107422 

-27.869299 

3.70 (3.58-3.82) 

10.46(10.24-10.47) 

101 (71-101) 

0.28 (0.10-0.28) 

222(45-240) 

z4 ERS 11888 

53.069328 

-27.714815 

3.70 (3.62-3.78) 

9.55 (9.29-9.76) 

49(19^9) 

0.00 (0.00-0.18) 

14 (14-77) 

z4 GSD 11269 

53.031239 

-27.785215 

3.71 (3.61-3.79) 

10.06(10.02-10.16) 

49 (40-101) 

0.10(0.10-0.24) 

45 (43-176) 

z4 GSD 27735 

53.138859 

-27.835371 

3.72 (3.65-3.79) 

9.56 (9.36-9.63) 

30(19-57) 

0.14(0.02-0.18) 

54 (17-78) 

z4 GSW 2898 

53.144775 

-27.871527 

3.78 (3.68-3.87) 

10.18(9.88-10.20) 

39 (30-49) 

0.20 (0.06-0.26) 

131 (35-200) 

z4 GSD 31543 

53.066261 

-27.849056 

3.81 (3.75-3.98) 

10 . 12 ( 10 . 00 - 10 . 20 ) 

49(30-57) 

0 . 12 ( 0 . 10 - 0 . 20 ) 

52(41-117) 

z4 GND 38889 

189.181412 

62.189281 

3.81 (3.75^.05) 

10.40(10.38-10.48) 

57 (49-71) 

0.14(0.14-0.22) 

68(65-163) 

z4 ERS 4079 

53.110340 

-27.689985 

3.81 (3.74-3.99) 

10.09 (9.87-10.17) 

19(10-57) 

0.32 (0.06-0.34) 

325 (31-373) 

z4 GSW 7015 

53.189873 

-27.892590 

3.83 (3.75-3.93) 

10.15 (9.77-10.15) 

49(30-80) 

0.16(0.06-0.22) 

84 (32-152) 

z4 GSW 6936 

53.073589 

-27.892235 

3.83 (3.76-3.92) 

10.47 (10.40-10.49) 

286 (202-718) 

0.00 ( 0 . 00 - 0 . 12 ) 

22 (22-81) 

z4 GSW 9851 

53.181850 

-27.906641 

3.85 (3.79-3.95) 

9.97 (9.64-9.94) 

49 (30-49) 

0.12 ( 0 . 00 - 0 . 20 ) 

56(18-114) 

z4 GNW 12987 

189.040390 

62.186352 

3.92 (3.84^.00) 

9.66 (9.49-9.69) 

57 (40-71) 

0.02 ( 0 . 00 - 0 . 02 ) 

19 (15-23) 

z4 ERS 5026 

53.133690 

-27.693453 

3.93 (3.82^.05) 

9.67 (9.63-9.90) 

39 (30-101) 

0.06 ( 0 . 00 - 0 . 12 ) 

27 (16^9) 

z4 GSD 34857 

53.076183 

-27.866360 

3.98 (3.89^.11) 

10.32(10.23-10.37) 

286(40-286) 

0.16 (0.14-0.16) 

123(96-128) 

z4 GSD 23593 

53.232452 

-27.822868 

3.99 (3.89^.14) 

9.48 (9.33-9.92) 

19(19-80) 

0.18(0.14-0.22) 

72 (52-101) 

z4 GSD 905 

53.168266 

-27.741940 

4.01 (3.92^.11) 

9.78 (9.44-9.76) 

90 (19-90) 

0.00 ( 0 . 00 - 0 . 20 ) 

13(13-88) 

z4 GNW 25081 

189.331039 

62.290836 

4.01 (3.90^.10) 

9.39 (9.34-9.43) 

10 ( 10 - 10 ) 

0.26 (0.24-0.26) 

140(122-144) 

z4 ERS 3543 

53.144386 

-27.687588 

4.02(3.87^.16) 

10.01 (9.96-10.26) 

39 (40^04) 

0.16(0.06-0.16) 

78(30-83) 

z4 GND 40720 

189.179291 

62.182003 

4.02 (3.84^.20) 

9.70 (9.26-9.74) 

30 (19-80) 

0.20 ( 0 . 02 - 0 . 20 ) 

75 (12-83) 

z4 GNW 18460 

189.286835 

62.367325 

4.05(3.92^.16) 

10.28 (10.22-10.38) 

286(202-806) 

0.00 ( 0 . 00 - 0 . 10 ) 

20 (19-53) 

z4 ERS 22264 

53.075272 

-27.755194 

4.06(3.96^.17) 

10.12(10.01-10.24) 

39 (30-508) 

0.16(0.10-0.16) 

93 (62-99) 

z4 GND 27047 

189.356354 

62.227554 

4.10(3.92^.21) 

10.47 (10.26-10.54) 

202 (49-904) 

0 . 12 ( 0 . 12 - 0 . 20 ) 

52 (50-98) 

z4 GNW 21799 

189.329178 

62.331532 

4.10(4.03^.17) 

9.82 (9.69-9.96) 

57 (40-80) 

0.00 (0.00-0.04) 

24 (23-35) 

z4 GNW 23907 

189.467804 

62.297764 

4.10(3.98^.21) 

9.94(9.82-10.31) 

80(19-101) 

0.20(0.14-0.24) 

122 (79-164) 
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TABLE At 

Stellar Populations of Bright Galaxies at z = 4 


z4 GND 7728 

189.270905 

62.291943 

4.10 

(4.01^.20) 

z4 GSD 15152 

53.027557 

-27.796583 

4.12 

(4.05^.22) 

z4 GNW 17778 

189.311584 

62.382095 

4.14 

(4.00^.29) 

z4 ERS 11615 

53.107601 

-27.713976 

4.14 

(4.01-4.26) 

z4 GND 23790 

189.233093 

62.236786 

4.15 

(4.06^.27) 

z4 GSD 20508 

53.192692 

-27.813051 

4.17 

(4.06^.26) 

z4 GNW 18340 

189.299255 

62.370079 

4.19 

(4.03^.30) 

z4 GNW 15232 

189.014038 

62.200378 

4.22 

(4.11-4.31) 

z4 ERS 14762 

53.016903 

-27.723013 

4.22 

(4.11^.32) 

z4 GNW 11056 

189.013092 

62.173153 

4.23 

(4.08^.31) 

z4 GSW 66 

53.230198 

-27.839573 

4.24 

(4.17^.33) 

z4 GSD 9138 

53.215435 

-27.778782 

4.24 

(4.10^.37) 

z4 GNW 3896 

189.145218 

62.129818 

4.24 

(4.12^.39) 

z4 GNW 8115 

189.145187 

62.154770 

4.25 

(4.09^.34) 

z4 GNW 26546 

189.445068 

62.281898 

4.26 

(4.02^.38) 

z4 GND 9704 

189.201706 

62.278599 

4.27 

(4.19^.36) 

z4 GND 21065 

189.290955 

62.244240 

4.29 

(4.13^.41) 

z4 GND 31952 

189.119812 

62.212612 

4.29 

(4.19^.37) 

z4 GNW 24624 

189.388474 

62.294460 

4.30 

(4.17^.43) 

z4 GNW 11698 

188.988419 

62.177494 

4.30 

(4.19^.39) 

z4 GNW 460 

189.120911 

62.101513 

4.30 

(4.18^.40) 

z4 GNW 2957 

189.185181 

62.122967 

4.30 

(4.17^.42) 

z4 GNW 31531 

189.295654 

62.349960 

4.31 

(4.15^.43) 

z4 GNW 7377 

189.108002 

62.151005 

4.31 

(4.25^.37) 

z4 GSD 21252 

53.143112 

-27.815502 

4.33 

(4.24^.39) 

z4 GND 30347 

189.324326 

62.217796 

4.34 

(4.24^.41) 

z4 GNW 24183 

189.466507 

62.296875 

4.35 

(4.21^.43) 

z4 GNW 21312 

189.366592 

62.330898 

4.38 

(4.25^.46) 

z4 GND 27301 

189.093216 

62.226814 

4.39 

(4.29^.47) 

z4 GSD 16522 

53.076065 

-27.800694 

4.39 

(4.32^.47) 

z4 GNW 21830 

189.366791 

62.331589 

4.40 

(4.30^.48) 

z4 GND 7158 

189.252747 

62.298897 

4.41 

(4.31^.50) 

z4 GNW 32075 

189.304321 

62.353088 

4.41 

(4.31^.48) 

z4 GNW 7816 

189.049530 

62.152920 

4.42 

(4.33^.49) 

z4 GSW 5453 

53.199787 

-27.884937 

4.42 

(4.32^.49) 

z4 GND 39360 

189.242157 

62.187439 

4.44 

(4.32^.54) 

z4 GSW 1319 

53.201015 

-27.860250 

4.44 

(4.37-4.51) 

z4 GNW 8416 

189.195557 

62.156864 

4.45 

(4.31^.54) 

z4 GND 32575 

189.173126 

62.210808 

4.45 

(4.36^.53) 

z4 GNW 13552 

189.061249 

62.189472 

4.45 

(4.34^.54) 

z4 GNW 8070 

189.107422 

62.154560 

4.45 

(4.37^.53) 

z4 ERS 5818 

53.069077 

-27.696518 

4.47 

(4.38^.56) 

z4 GNW 20572 

189.313507 

62.320702 

4.47 

(4.37-4.56) 

z4 GNW 18301 

189.335205 

62.370796 

4.47 

(4.39^.54) 

z4 GNW 18575 

189.297821 

62.365429 

4.49 

(4.42^.55) 

z4 GNW 2261 

189.192337 

62.117081 

4.50 

(4.42-4.57) 

z4 GND 40010 

189.335739 

62.184937 

4.51 

(4.44^.56) 

z4 GNW 7213 

189.108322 

62.150093 

4.51 

(0.69^.57) 

z4 GNW 7206 

189.108185 

62.149807 

4.54 

(4.46^.61) 

z4 GSW 8512 

53.177471 

-27.900093 

4.54 

(4.48^.60) 

z4 GND 25942 

189.147903 

62.230583 

4.54 

(4.44^.63) 

z4 GSD 36028 

53.079254 

-27.877260 

4.54 

(4.46^.59) 

z4 GND 14271 

189.305084 

62.263287 

4.56 

(4.44^.65) 

z4 GNW 18613 

189.276077 

62.364826 

4.57 

(4.47^.64) 


10.23 (9.78-10.24) 
10.41 (10.23-10.44) 
9.76 (9.29-9.82) 

9.52 (9.61-9.96) 
10.50(10.49-10.60) 

9.65 (9.61-9.92) 
11.29(11.17-11.33) 
10.18(9.80-10.18) 

9.40 (9.32-9.84) 
9.64 (9.61-9.78) 
9.76 (9.76-10.22) 
10.07 (9.82-10.08) 

9.80 (9.77-10.12) 
9.62 (9.54-9.79) 

10.36(10.35-10.59) 

9.85 (9.85-10.23) 
10.47 (10.17-10.53) 

9.56 (9.50-9.90) 
9.34 (9.35-9.62) 
9.62 (9.42-9.66) 

9.41 (9.36-9.80) 
9.90(9.89-10.12) 

10.32(10.10-10.46) 
10.35 (10.24-10.35) 
10.05 (10.03-10.47) 
9.71 (9.54-9.82) 
10.09 (9.95-10.22) 

9.66 (9.72-10.34) 
10.00 (9.90-10.14) 

9.67 (9.52-10.14) 

9.95 (9.95-10.26) 
10.27 (9.94-10.30) 
10.46(10.45-10.77) 

9.70 (9.71-10.15) 
9.90 (9.77-10.03) 

9.86 (9.79-10.15) 
9.84 (9.60-9.87) 

10.51 (10.39-10.58) 
9.75 (9.37-9.72) 
9.92 (9.94-10.38) 

9.81 (9.71-9.88) 
9.43 (9.37-9.64) 
9.83 (9.52-9.86) 
9.97 (9.88-10.31) 

10.70(10.37-10.74) 

9.96 (9.62-9.97) 
10.16(9.98-10.20) 

9.89 (6.51-9.87) 

9.97 (9.85-10.12) 

9.53 (9.54-10.02) 
10.09 (9.89-10.23) 
10.21 (9.96-10.22) 
10.79(10.45-10.79) 
10.04(10.06-10.58) 


101 (19-101) 
101(49-101) 
30(19^9) 
19(30-101) 
101 (57^04) 
30 (30-202) 
1015(904-1015) 
71 (30-80) 

19 (19-80) 

39 (30-57) 
19(19-90) 
101 (80-101) 
10 (10-40) 

49(40-80) 
101 (57-286) 
30(30-90) 
202(80-286) 
30 (30-57) 

39 (30-57) 

30 (19^9) 
30(19-71) 

30 (30-49) 
404(49-286) 
80(30-80) 

57 (30-101) 
1015 (40-404) 
570 (71-570) 
10(10-30) 
904(40-718) 
19(10-90) 

30 (30-286) 
101 (40-101) 
30 (30-202) 
19(19-71) 
286(80-570) 
71(30-90) 

39 (30-90) 
286(49-286) 
39 (30-49) 
19(19-71) 

30 (30-49) 

39 (30-49) 

39(30-90) 

30 (30-90) 
30(19-57) 

39(30-80) 

57(30-90) 
30(19-90) 
19(19-30) 
19(19^0) 
19(19-57) 
1015 (30-202) 
39 (40-80) 

19 (19-57) 


0.20 (0.16-0.22) 
0.20 (0.20-0.24) 
0.18 (0.10-0.20) 
0.18 (0.02-0.18) 
0.18(0.08-0.18) 
0.14(0.04-0.16) 
0.18(0.18-0.24) 
0.22 (0.20-0.26) 
0.16(0.08-0.18) 
0.04 (0.00-0.08) 
0.18(0.12-0.18) 
0.18 (0.00-0.20) 
0.34 (0.32-0.34) 
0.02 (0.00-0.06) 
0.28 (0.00-0.30) 
0.14(0.04-0.16) 
0.30 (0.18-0.30) 
0.12 (0.06-0.14) 
0.00 (0.00-0.10) 
0.16(0.02-0.20) 
0.08 (0.00-0.16) 
0.16(0.12-0.18) 
0.24(0.16-0.30) 
0.00 (0.00-0.12) 
0.18 (0.16-0.22) 
0.08 (0.06-0.12) 
0.14(0.10-0.18) 
0.26 (0.22-0.32) 
0.14(0.10-0.18) 
0.22 (0.22-0.26) 
0.20 (0.14-0.20) 
0.26 (0.22-0.28) 
0.26 (0.20-0.26) 
0.18(0.12-0.20) 
0.16(0.14-0.20) 
0.18(0.16-0.24) 
0.14(0.00-0.18) 
0.24 (0.20-0.30) 
0.14(0.06-0.18) 
0.24 (0.20-0.26) 
0.16(0.04-0.18) 
0.02 (0.00-0.16) 
0.16(0.04-0.16) 
0.18 (0.16-0.22) 
0.28 (0.24-0.30) 
0.14(0.02-0.18) 
0.06 (0.00-0.16) 
0.20 (0.10-0.66) 
0.26 (0.22-0.28) 
0.18(0.18-0.26) 
0.30 (0.24-0.32) 
0.10(0.06-0.12) 
0.36 (0.30-0.36) 
0.22 (0.20-0.26) 


133(91-163) 
171 (178-266) 
85(41-99) 
79(18-83) 
211 (83-221) 

59 (25-76) 
146(142-226) 
178(139-242) 

57 (28-78) 
26 (16-38) 
138(80-143) 
92(18-102) 
346 (297-359) 
19(15-31) 
151 (12-228) 
95(38-117) 
264(81-240) 
45 (27-58) 
12 ( 11 ^ 1 ) 
61(14-82) 
27(14-64) 
84 (59-95) 
122 (54-208) 
57(56-179) 
181(149-266) 
26 (21-39) 
68(42-107) 
210(154-394) 
50(37-74) 
160 (154-224) 
118(65-120) 
146 (90-167) 
384 (229-397) 
100(64-121) 
52(42-75) 
65 (54-118) 

60 (16-89) 
136(91-259) 

49 (21-72) 
164(109-184) 
96 (29-105) 
16(13-65) 
58(17-56) 
116(98-194) 
739(458-894) 
79 (23-115) 
62 (35-158) 
116(0-142) 
245(149-279) 
67 (67-152) 
325(180-360) 
83 (56-104) 
537(295-541) 
219(184-335) 
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TABLE A2 

Stellar Populations of Bright Galaxies at z = 5 



ID 

Right Ascension 
(J2000) 

Declination 

(J2000) 

Redshift 

log M* 

(Mq) 

Age 

(Myr) 

E(B-V) 

SER 

(Mo yr‘) 

z5 

GSW 8762 

53.208008 

-27.901289 

4.51 (4.44^.57) 

9.72 (9.50-9.78) 

49 (30-71) 

0.06 (0.00-0.08) 

31(16-38) 

z5 

GND 36639 

189.186142 

62.197327 

4.56 (4.45^.68) 

10.25 (10.11-10.52) 

80(19-101) 

0.28 (0.26-0.32) 

223 (203-330) 

z5 

GNW 17976 

189.300125 

62.377483 

4.57 (4.42-4.74) 

10.19(10.12-10.37) 

286(80-286) 

0.22(0.10-0.26) 

102 (31-164) 

z5 

GSD 9044 

53.091724 

-27.778580 

4.58 (4.48^.67) 

10.26 (10.02-10.25) 

202 (40-508) 

0.06(0.06-0.18) 

32 (33-110) 

z5 

GSW 6918 

53.234589 

-27.892109 

4.66 (4.58^.75) 

9.52 (9.38-9.59) 

57 (49-57) 

0.00 (0.00-0.02) 

14(13-18) 

z5 

GND 38041 

189.299210 

62.192570 

4.67 (4.60^.74) 

9.36 (9.31-9.89) 

10 (10-19) 

0.22 (0.20-0.28) 

106 (94-198) 

z5 

GND 16758 

189.176086 

62.256329 

4.71 (4.62^.84) 

10.03 (9.96-10.11) 

19 (19-30) 

0.26 (0.26-0.28) 

277 (272-336) 

z5 

GSD 33149 

53.070778 

-27.856453 

4.71 (4.62^.81) 

10.02 (9.63-10.30) 

101 (10-508) 

0.20(0.10-0.24) 

116(49-172) 

z5 

GSD 13326 

53.095345 

-27.790989 

4.71 (4.64^.79) 

10.19(9.98-10.26) 

30 (19-80) 

0.26(0.16-0.26) 

228(86-238) 

z5 

i ERS 3475 

53.070839 

-27.687143 

4.72(4.66^.81) 

9.95 (9.84-10.39) 

19(19^0) 

0.26 (0.26-0.32) 

205 (211-369) 

z5 

GND 38212 

189.273590 

62.192028 

4.72 (4.65-5.06) 

9.84 (9.82-10.19) 

30 (30-71) 

0.16(0.08-0.20) 

72 (43-136) 

z5 

GSW 1565 

53.238213 

-27.862486 

4.75 (4.70^.79) 

9.96 (9.51-9.96) 

30 (19-30) 

0.18(0.12-0.20) 

135(73-158) 

z5 

GND 12253 

189.294418 

62.269447 

4.77 (4.70^.90) 

10.06 (9.87-10.19) 

570(40-718) 

0.14(0.10-0.20) 

55 (39-104) 

z5 

ERS 12604 

53.021912 

-27.716784 

4.81 (4.76^.89) 

9.60 (9.60-9.84) 

39(40-1015) 

0.08 (0.02-0.10) 

30 (17-34) 

z5 

GNW 25539 

189.489624 

62.288536 

4.81 (4.76^.87) 

10.34 (10.03-10.43) 

19 (19^0) 

0.26(0.18-0.26) 

570(266-567) 

z5 

i ERS 2517 

53.119019 

-27.682158 

4.82 (4.76^.93) 

9.70 (9.41-9.77) 

30 (30-80) 

0.18(0.00-0.16) 

73(14-68) 

z5 

GND 17343 

189.091400 

62.254662 

4.83 (4.75^.95) 

10.16(10.05-10.19) 

30 (19-71) 

0.24 (0.06-0.24) 

214(45-230) 

z5 

i ERS 6044 

53.048820 

-27.697111 

4.85 (4.78^.94) 

9.85 (9.78-9.87) 

90 (30-90) 

0.00 (0.00-0.20) 

15 (15-94) 

z5 

GSW 6966 

53.245884 

-27.892273 

4.94 (0.83^.98) 

10.29 (10.20-10.37) 

19 (10-30) 

0.34 (0.28-0.36) 

506 (265-566) 

z5 

GNW 6112 

189.064835 

62.143963 

4.95 (4.90-5.06) 

9.85 (9.47-10.02) 

19 (10-71) 

0.26 (0.06-0.28) 

183(31-221) 

z5 

GSD 10352 

53.021172 

-27.782366 

4.95 (4.91^.99) 

10.43 (10.36-10.80) 

57 (19-90) 

0.22 (0.08-0.24) 

474(137-571) 

z5 

GNW 4779 

189.203064 

62.136204 

4.96 (4.88-5.15) 

9.86(10.07-10.31) 

19 (30-286) 

0.28 (0.04-0.32) 

174 (20-244) 

z5 

GND 15230 

189.205124 

62.260712 

5.02 (4.96-5.08) 

9.41 (9.41-9.69) 

39 (30-49) 

0 .00(0.00-0.12) 

15(15-59) 

z5 

GNW 16101 

189.000168 

62.207241 

5.02 (4.78-5.06) 

9.50 (9.20-9.52) 

49 (30-49) 

0.02 (0.00-0.04) 

18(12-20) 

z5 

GNW 10657 

189.130219 

62.170780 

5.03 (4.80-5.15) 

9.35 (9.19-9.44) 

49 (30-57) 

0.00 (0.00-0.06) 

13 (11-27) 

z5 

GNW 19973 

189.371780 

62.324554 

5.03 (4.80-5.12) 

10.10(10.15-10.59) 

10 (10-202) 

0.38 (0.22-0.34) 

681 (140^73) 

z5 

i ERS 2314 

53.088764 

-27.680889 

5.05 (5.01-5.20) 

9.55 (9.43-9.60) 

30 (19^9) 

0.14(0.00-0.20) 

52 (16-86) 

z5 

PARI 1735 

53.246998 

-27.686445 

5.07 (4.99-5.18) 

10.06 (9.74-10.09) 

49 (30-71) 

0.16(0.08-0.18) 

69 (32-79) 

z5 

GNW 13254 

189.039703 

62.187725 

5.24(5.15-5.32) 

9.88 (9.74-9.93) 

71(49-71) 

0.00 (0.00-0.04) 

22 (21-33) 

z5 

i ERS 9511 

53.155231 

-27.700718 

5.26 (5.05-5.40) 

9.90 (9.49-9.99) 

49 (19-90) 

0.12 (0.00-0.20) 

47(14-98) 

z5 

GNW 25408 

189.457428 

62.289562 

5.27 (5.11-5.40) 

9.62 (9.35-9.66) 

71 (40-80) 

0.00 (0.00-0.02) 

12 (11-15) 

z5 

PARI 1385 

53.273891 

-27.683685 

5.29 (5.20-5.36) 

9.66 (9.67-9.99) 

30 (30-71) 

0.16(0.12-0.22) 

60 (42-113) 

z5 

GSD 17901 

53.200577 

-27.804909 

5.33 (5.23-5.43) 

9.91 (9.57-9.97) 

57 (30-90) 

0.08(0.00-0.12) 

34 (17^8) 

z5 

GND 34380 

189.249146 

62.205204 

5.35 (5.24-5.45) 

10.42 (10.35-10.61) 

49 (40-508) 

0.18(0.14-0.24) 

136(85-234) 

z5 

GND 33094 

189.266449 

62.209164 

5.36 (5.26-5.43) 

9.82 (9.54-9.79) 

49 (40-80) 

0 .10(0.02-0.10) 

39 (16-36) 

z5 

GNW 27194 

189.523804 

62.278839 

5.40 (5.30-5.50) 

10.14(10.01-10.30) 

57 (40-202) 

0.08 (0.04-0.16) 

50(33-105) 

z5 

GND 39570 

189.175720 

62.186714 

5.45 (5.33-5.55) 

9.74 (9.80-10.22) 

19 (30-49) 

0 .22(0.12-0.22) 

126(49-140) 

z5 

GND 10054 

189.169846 

62.277298 

5.46 (5.37-5.53) 

9.92 (9.62-9.91) 

57 (40-90) 

0.08 (0.00-0.14) 

35 (16-58) 

z5 

GND 35096 

189.283508 

62.203049 

5.48 (5.35-5.55) 

10.32 (9.81-10.38) 

202 (19-101) 

0.08 (0.04-0.22) 

43(31-167) 

z5 

GSD 4579 

53.170231 

-27.762848 

5.48 (5.42-5.54) 

9.67 (9.49-9.76) 

30 (19^9) 

0.12(0.08-0.14) 

49(33-60) 

z5 

GNW 3960 

189.055939 

62.129990 

5.52 (5.38-5.67) 

9.88 (9.59-9.88) 

71(30-90) 

0.02 (0.00-0.08) 

23 (16^2) 

z5 

GND 18617 

189.353439 

62.250774 

5.53 (5.46-5.63) 

9.68 (9.56-9.91) 

19 (19-30) 

0.16(0.12-0.18) 

95 (67-122) 

z5 

GND 24948 

189.320312 

62.233444 

5.54 (5.47-5.63) 

9.53 (9.54-10.10) 

10 (10-30) 

0.26 (0.26-0.32) 

179(178-331) 

z5 

GNW 28218 

189.503250 

62.273884 

5.55 (5.48-5.62) 

10.47 (10.16-10.47) 

30 (19^9) 

0.26(0.14-0.26) 

435 (125^31) 

z5 

GSD 4436 

53.170811 

-27.762228 

5.58 (5.49-5.65) 

9.53 (9.53-9.78) 

30 (19^9) 

0.10(0.02-0.18) 

35(17-84) 

z5_ 

GNW 29490 

189.312943 

62.344555 

5.69 (5.55-5.84) 

9.58 (9.58-9.89) 

10 (10-19) 

0.28 (0.28-0.30) 

200 (182-254) 
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Increasing Stellar Baryon Fractions at z > 4 

TABLE A3 

Stellar Populations of Bright Galaxies at z = 6 


ID 

Right Ascension 
(J2000) 

Declination 

(J2000) 

Redshift 

log M* 

(Mo) 

Age 

(Myr) 

E(B-V) 

SER 

(M 0 yr-‘) 

z 6 GSD 29074 

53.156788 

-27.839560 

5.56 

(5.46-5.66) 

9.71 (9.54-10.02) 

30(19-90) 

0.16(0.00-0.18) 

67(16-85) 

z 6 GSD 17919 

53.074165 

-27.804928 

5.65 

(5.53-5.75) 

9.72 (9.55-9.85) 

80 (57-202) 

0.00 (0.00-0.04) 

13 (12-20) 

z 6 GNW 10970 

189.075775 

62.172729 

5.67 

(5.57-5.75) 

10.52 (10.37-10.56) 

404 (40-570) 

0.20 (0.18-0.24) 

170(146-243) 

z 6 GND 16399 

189.234833 

62.257507 

5.67 

(5.58-5.76) 

9.66 (9.44-9.75) 

39 (19-57) 

0.08 (0.02-0.18) 

27 (16-87) 

z 6 GND 16819 

189.328232 

62.256317 

5.69 

(5.61-5.76) 

10.09 (9.82-10.20) 

19 (19-30) 

0.28 (0.20-0.26) 

319(130-302) 

z 6 ERS 4104 

53.066242 

-27.689983 

5.72 

(5.61-5.84) 

9.57 (9.48-9.67) 

39 (30^9) 

0.02 (0.00-0.04) 

21 (17-27) 

z 6 GSW 12831 

53.106689 

-27.930193 

5.79 

(5.72-5.84) 

10.05 (9.96-10.24) 

49(40-101) 

0.06 ( 0 . 00 - 0 . 10 ) 

57 (34-83) 

z 6 GSW 6659 

53.151745 

-27.890762 

5.80 

(5.68-5.91) 

9.65 (9.39-9.82) 

39(19-71) 

0.08 (0.00-0.14) 

31(14-52) 

z 6 GSD 27934 

53.101673 

-27.836084 

5.80 

(5.70-5.90) 

10.61 (10.59-10.67) 

19 (19-30) 

0.38 (0.36-0.40) 

1066 (863-1186) 

z 6 MAIN 5871 

53.166721 

-27.804167 

5.80 

(5.76-5.83) 

9.90 (9.64-9.91) 

57 (40-57) 

0.04 (0.00-0.04) 

34 (21-34) 

z 6 GND 28043 

189.418427 

62.224796 

5.85 

(5.75-5.94) 

10.00(9.69-10.14) 

49 (19-80) 

0.10(0.02-0.18) 

46 (25-104) 

z 6 GNW 22717 

189.416199 

62.333141 

5.88 

(5.66-6.04) 

10.17 (9.74-10.27) 

806 (30-806) 

0.20 ( 0 . 02 - 0 . 22 ) 

94(17-117) 

z 6 GSD 23051 

53.225368 

-27.821125 

5.96 

(5.86-6.05) 

10.07 (9.91-10.26) 

39 (30-202) 

0.14(0.02-0.18) 

89 (30-122) 

z 6 GNW 23437 

189.388000 

62.301167 

6.02 

(5.85-6.18) 

9.84 (9.60-9.87) 

101 (57-101) 

0.00 (0.00-0.08) 

12(11-27) 

z 6 PARI 1068 

53.234592 

-27.680861 

6.13 

(6.05-6.23) 

10.10(9.90-10.16) 

39 (40-80) 

0.18 (0.06-0.18) 

109(34-103) 

z 6 GNW 22555 

189.301788 

62.307808 

6.15 

(6.06-6.25) 

9.63 (9.41-9.65) 

57(40-57) 

0.02 ( 0 . 00 - 0 . 02 ) 

18(12-19) 

z 6 ERS 7413 

53.158138 

-27.702112 

6.25 

(6.12-6.39) 

9.91 (9.66-10.07) 

30 (19^9) 

0.10(0.00-0.14) 

85(38-127) 

z 6 GSD 21289 

53.076241 

-27.815453 

6.30 

(6.20-6.39) 

9.63 (9.54-9.88) 

39 (30-202) 

0.10(0.00-0.16) 

32(13-63) 

z 6 GSD 233 

53.188583 

-27.733210 

6.39 

(6.29-6.48) 

9.88 (9.42-9.93) 

30 (19^9) 

0.20 (0.04-0.20) 

112(23-115) 


TABLE A4 

Stellar Populations of Bright Galaxies at z = 7 


ID 

Right Ascension 
(J2000) 

Declination 

(J2000) 

Redshift 

log Mt 
(Mq) 

Age 

(Myr) 

E(B-V) 

SER 

(Mq yr-‘) 

z7 

GSD 25074 

53.233047 

-27.827383 

6.66 (6.47-6.79) 

10.09(9.87-10.15) 

404 (49-508) 

0.18 (0.00-0.14) 

82(15-55) 

z7 

ERS 12574 

53.094410 

-27.716846 

6.71 (6.58-6.85) 

9.83 (9.63-10.06) 

49 (30-202) 

0.08 (0.00-0.18) 

35 (16-90) 

z7 

GNW 24443 

189.356888 

62.295319 

6.72 (6.65-6.81) 

10.02 (9.80-10.05) 

30 (30-71) 

0.18 (0.00-0.18) 

155(29-160) 

z7 

GNW 24671 

189.361710 

62.294373 

6.72 (6.27-7.07) 

9.15 (8.53-9.25) 

19 (10-40) 

0 . 10 ( 0 . 00 - 0 . 10 ) 

37 (11-33) 

z7 

GSD 21368 

53.154922 

-27.815744 

6.75 (6.62-6.97) 

10.14(9.89-10.34) 

404 (30-286) 

0.20 (0.00-0.26) 

82 (14-142) 

z7 

GNW 17001 

189.032486 

62.216415 

6.76 (6.55-6.92) 

10.74 (10.70-10.82) 

404 (286-570) 

0.06 ( 0 . 02 - 0 . 12 ) 

22 (14-37) 

z7 

GSD 21172 

53.155342 

-27.815178 

6.84 (6.76-6.96) 

10.19(9.95-10.20) 

30 (19-80) 

0.24 (0.00-0.30) 

231(24-396) 

z7 

GSD 10175 

53.210335 

-27.782211 

7.05 (6.78-7.32) 

9.70 (9.23-9.81) 

19 (10-71) 

0.22 (0.00-0.24) 

129(15-141) 

z7 

GND 18181 

189.082687 

62.252476 

7.09 (6.94-7.28) 

9.36 (9.28-9.96) 

10 (10-49) 

0.22 (0.12-0.24) 

123(47-142) 

z7 

GND 11402 

189.186172 

62.270863 

7.16(7.00-7.34) 

9.69 (9.22-9.76) 

30 (10-49) 

0.10(0.02-0.18) 

51(21-114) 

z7 

GNW 4703 

189.094528 

62.135540 

7.17(7.04-7.37) 

9.55 (9.06-9.68) 

39 (10-49) 

0.02 ( 0 . 00 - 0 . 10 ) 

27 (22-61) 

z7 

GNW 19939 

189.273392 

62.324783 

7.25 (7.03-7.75) 

8.79 (8.64-9.29) 

10 (10-30) 

0 . 10 ( 0 . 00 - 0 . 12 ) 

34(14^2) 

z7 

GND 42912 

189.157883 

62.302372 

7.49 (7.33-7.70) 

9.51 (9.41-9.65) 

10 ( 10 - 10 ) 

0.24 (0.20-0.28) 

178(133-224) 

z7 

PAR2 3098 

53.281712 

-27.867699 

7.66 (7.44-7.82) 

9.82 (9.56-9.90) 

30 (19^9) 

0.14(0.02-0.18) 

69(23-111) 










